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Abstract 

We investigate the mass spectrum of the pentaquark baryon (ududs) 
in quenched lattice QCD with exact chiral symmetry. Using 3 different 
interpolating operators for G, we measure the 3x3 correlation matrix 
and obtain the eigenvalues A^(t) with ± parity, for 100 gauge configu- 
rations generated with Wilson gauge action at (3 = 6.1 on the 20 3 x 40 
lattice. For the lowest-lying J p = l/2~ state, its mass is almost identi- 
cal to that of the KN s-wave, while for the lowest-lying J p = l/2 + state, 
its mass is smaller than that of the KN p-wave, especially for the regime 
m u < m s . By chiral extrapolation (linear in m^) to = 135 MeV, 
we obtain the masses of the lowest-lying states: m(l/2~) = 1424(57) 
MeV and m(l/2 + ) = 1562(121) MeV, in agreement with the masses of 
rriK + rriN ~ 1430 MeV and + (154O) respectively. 
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1 Introduction 



The recent experimental observation of the exotic baryon + (154O) (with the 
quantum numbers of K + n) by LEPS collaboration P at Spring-8 and the sub- 
sequent confirmation [21 El El El El El E] from some experimental groups has 
become one of the most interesting topics in hadron physics. The remarkable 
features of O + (1540) are its strangeness S = +1, and its exceptionally narrow 
decay width (< 15 MeV) even though it is ~ 100 MeV above the KN thresh- 
old. Its strangeness S = +1 immediately implies that it cannot be an ordinary 
baryon composed of three quarks. Its minimal quark content is ududs. Nev- 
ertheless, there are quite a number of experiments jH] which so far have not 
observed G + (1540) or any pentaquarks. This casts some doubts about the 
existence of 6 + (1540). 

Historically, the experimental search for G + (1540) was motivated by the 
predictions of the chiral-soliton model^U], an outgrowth of the Skyrme modelfTT] 
Even though the chiral solition model seems to provide very close predictions 
for the mass and the width of O + (1540), obviously, it cannot reproduce all 
aspects of QCD, the fundamental theory of strong interactions. Now the cen- 
tral theoretical question is whether the spectrum of QCD possesses + with the 
correct quantum numbers, mass, and decay width. 

At present, the most viable approach to solve QCD nonperturbatively from 
the first principles is lattice QCD. Explicitly, one needs to construct an inter- 
polating operator which has a significant overlap with the pentaquark baryon 
states. Then one computes the time-correlation function of this interpolating 
operator, and from which to extract the masses of its even and odd parity states 
respectively. However, any (ududs) operator must couple to hadronic states 
with the same quantum numbers (e.g., KN scattering states). It is neces- 
sary to disentangle the lowest-lying pentaquark states from the KN scattering 
states, as well as the excited pentaquark states. 

To this end, we adopt the so-called variational method [121 EH] , and use 
three different interpolating operators for Q(ududs), to compute their 3x3 
correlation matrix, and from its eigenvalues to extract the masses of the even 
and odd parity states. This is the first lattice QCD study of G + with 3x3 
correlation matrix. These three interpolating operators (with 1 = 0) are: 

[u T C^ 5 d] xc {s x/ 3 e (7 5 ) ( g r ,u :!;r , e (7 5 d) a . Q , c - 8^(75)^^(7511)3.0,,,} (1) 
[u C7 5 d] xc {s x/3e (75) 

/3rj^-xric (75 d) xae &x/3e (T5) f3rjd-xr)c{ r )b^) xae }(2) 

e cde [u T C-f 5 d} xc [u T Cd} xd (Cs T ) xae (3) 

where u, d and s denote the quark fields; e c d e is the completely antisymmetric 
tensor; x, {c, d, e) and {a, j3, rj} denote the lattice site, color, and Dirac indices 
respectively; and C is the charge conjugation operator satisfying Oy^C -1 = 



(01) xa = 

(02) xa = 

(Os) XQ = 



1 



— 7J and {C^) T = — C75. Here the diquark operator is defined as 

[u T rd] xa = e a b c (u xab T 'apd-xpc — d xab T af3 u x i3 C ) (4) 

where r Q/ g = — Tp a . Thus the diquark transforms like a spin singlet (l s ), color 
anti-triplet (3 C ), and flavor anti-triplet (3/). For T = C75, it transforms as a 
scalar, while for r = C, it transforms like a pseudoscalar. Here 0\, O2, and 
O3 all transform like an even operator under parity. 

The operator 0\ is similar to the naive kaon ® nucleon operator which was 
used by Mathur et al. [Tl]. The operator O2 was first considered by Zhu ^3] 
and was adopted by Csikor et al. ^H] i n their lattice study. The difference 
between 0\ and O2 is that, in the latter case, the color index of the u(d) quark 
in the kaon is swapped with that of the d(u) quark in the nucleon such that 
the "kaon" and the "nucleon" do not appear as color singlets. The operator 
O3 is motivatied by the Jaffe- Wilzcek (diquark-diquark-antiquark) model [T7J , 
which was considered by Sugiyama et al. [TH] . Sasaki [TH], and was adopted 
by Chiu and Hsieh [20], and Ishii et al. [2T] . 

In the Jaffe- Wilzcek model, each pair of [ud] form a diquark which trans- 
forms like a spin singlet (l s ), color anti-triplet (3 C ), and flavor anti-triplet 
(3/). Then the pentaquark baryon 0([txd][ixd]s) emerges as the color singlet 
in (3 C x 3 C ) x 3 C = l c + 8 C + 8 C + 10 c , and a member (with S — +1 and 
I = 0) of the flavor anti-decuplet in 3/ x 3/ x 3/ = 1/ + 8/ + 8/ + 10/. Now, 
if one attempts to construct a local interpolating operator for [ud][ud]s, then 
these two identical diquarks must be chosen to transform differently (i.e., one 
scalar and one pseudoscalar), otherwise e a b c [ud] x b[ud] xc s xaa is identically zero 
since diquarks are bosons. Thus, when the orbital angular momentum of this 
scalar-pseudoscalar-antifermion system is zero (i.e., the lowest lying state), its 
parity is even rather than odd. Alternatively, if these two diquarks are located 
at two different sites, then both diquark operators can be chosen to be scalar, 
however, they must be antisymmetric in space, i.e., with odd integer orbital 
angular momentum. Thus the parity of lowest lying state of this scalar-scalar- 
antifermion system is even, as suggested in the original Jaffe- Wilzcek model. 
(Note that all correlated quark models e.g., Karliner-Lipkin model [21] and 
flavor-spin model (22]; advocate that the parity of O + (1540) is positive.) 

Evidently, the diquark operator plays an important role in constructing 
sources for pentaquark baryons as well as 3-quark baryons. The possibility 
of forming multiquark hadrons through diquark correlations was proposed by 
Jaffe in 1977 [2H]. Although the idea is essentially based on the color-spin 
interaction between the quarks (through one gluon exchange), its salient fea- 
tures seem to persist even at the hadronic distance scale where QCD is strongly 
coupled. Thus, it is interesting to see whether such multiquark hadrons (e.g., 
pentaquark baryons) do exist in the spectrum of QCD. 

In this paper, we use the optimal domain-wall fermion [27] to study the pen- 
taquark baryons. The salient features of optimal lattice domain-wall fermion 
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are: (i) The quark propagator as well as the effective 4D lattice Dirac operator 
for internal fermion loops have optimal chiral symmetry for any N s (number 
of sites in the 5-th dimension) and gauge background; (ii) The quark fields 
and hadron observables manifest the discrete symmetries of their counterparts 
in continuum; (iii) The quark action is ultralocal on the 5 dimensional lat- 
tice, thus the dynamical quark can be simulated with the standard hybrid 
monte carlo algorithm; (iv) The quark propagator in gauge background can 
be computed efficiently through the effective 4D lattice Dirac operator. 

Before we turn to our results, it is instructive to review the current status 
of quenched lattice QCD for + , as summarized in Table ^ At first sight, 
current lattice QCD results for + seems to disagree with each other. Note 
that the different claims listed in the last column of Table ^ already cover 
all possible outcomes. Obviously, some of the claims in Table [0 cannot be 
sustained for a long time, no matter what is the experimental outcome. Now 
if one compares the essential features among these exploratory lattice studies, 
one might understand what could be the causes for these different claims. In 
the following, we pinpoint the crucial features which may have direct impacts 
to these claims. 

So far, all lattice QCD simulations are performed at unphysically large 
m u , thus it is necessary to chirally extrapolate to physical m u (or equivalently 
m n = 135 MeV). Then a crucial question is how good are the data points 
used for chiral extrapolation, i.e., how many data points are obtained with 
m u < m s (an obviously physical condition ought to be satisfied), and what is 
the smallest m u (or equivalently m n ) in these data points. These two questions 
are answered in the columns with headings N(m u < m s ) and m!^ %n respectively. 
Obviously, if m™ m is too large, and/or the number N(m u < m s ) is too small, 
then the chiral extrapolation would tend to overestimate the masses, especially 
for the excited states. 

Another important question is whether the interpolating operator one uses 
has a significant overlap with the pentaquark state. If it has little overlap with 
the pentaquark state, then the signal might be too weak to be detected. As 
we will see below, 0±, O2 and O3 all have good overlap with the lowest-lying 
negative parity state for the entire range of m u . However, the lowest-lying 
negative parity state turns out to be nothing but the KN s-wave scattering 
state. On the other hand, for the positive parity channel, only O3 has a 
significant overlap with the lowest-lying positive parity state, in the regime 
m u < m s /2. As we will see below, this positive parity state is ruled out to 
be KN p-wave, KN* s-wave (where N* is the lowest negative parity state 
of nucleon), or KNrj' s-wave (where rf is the artifact due to the quenched 
approximation) . 

At the end of this paper, we will return to Tabled to discuss what could 
be the causes for the different claims in these exploratory lattice studies. 
The outline of this paper is as follows. In Section 2, we outline our com- 
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Ref. 


Operator 


Quark 


Gauge 


a 


L 


m min 


N(m u < m s ) 


m(l/2-), S/R 


m(l/2+), S/R 


Signal/Parity 




Oi + a0 2 


W 


W 


0.09 


1.8 


420 


3 


1539(50), R 


2710(79) 


Yes/- 


nsi 


o 3 


W 


W 


0.07 


2.2 


650 


1 


1840(80), R 


2940(130) 


Yes/- 


This work 


{Cl 5 2,3}3x3 


Odwf 


W 


0.09 


1.8 


440 


10 


1433(72), S 


1562(121), R 


Yes/+ 


nu 


o 1 


Ov 


Iw 


0.20 


2.4/3.2 


180 


13 


~ 1450, S 


~ 1650, S 


No 




o 3 


W 


W 


0.18 


2.1 


656 





1750(40), S 


2250(120) 


No 



Table 1: Summary of current lattice QCD results for G + . Here all masses are in units of MeV, the lattice spacing (a) and 
the box size (L) are in units of Fermi. The notations for the symbols are: W - Wilson fermion or Wilson gauge action; 
Ov - overlap fermion; Iw - Iwasaki gauge action; m™ m - the smallest pion mass; N(m u < m s ) - the number of data 
points satisfying the condition m u < m s ; m(l/2 ± ) - the mass of l/2 ± state via chiral extrapolation; S/R - Scattering 
state/Resonance. There are two lattice studies [221I2SI which have not listed here, due to lack of the information for some 
of the entries. 



putatation of quark propagators. In Section 3, we outline our determination 
of the lattice spacing a, and the strange quark bare mass m s . In Section 4, we 
present our results of the masses of the even and odd parity states extracted 
from the 3x3 correlation matrix of Oi, O2 and O3. These are the first lattice 
QCD results using 3x3 correlators for + . In Section 5, we investigate the 
KN scattering states with "disconnected" KN (i.e., without quark exchanges 
between K and N) correlation function, and use them to identify the KN scat- 
tering states in the spectrum of 3 x 3 correlation matrix. In Section 6, we 
discuss the current lattice QCD results for + , and conclude with some re- 
marks. In the Appendix, we include our results of the masses of the even and 
odd parity states extracted from the time-correlation functions of Oi, O2, and 
3 respectively. 



2 Computation of quark propagators 

Now it is straightforward to work out the baryon propagator (® X a®ys) in 
terms of quark propagators. In lattice QCD with exact chiral symmetry, quark 
propagator with bare mass m q is of the form (D c + m q )~ l [2E]> where D c is 
exactly chirally symmetric at finite lattice spacing. In the continuum limit, 
(D c +m 5 ) _1 reproduces [ r y iJi {d^+iA i _ l )+m q ]~ l . For optimal domain- wall fermion 
with N s + 2 sites in the fifth dimension, 

1 + j 5 S(H w ) 
D c = 2m 0z 



S(H U 



1-1 5 S(H U 

i-nf=iT s 



i + nfii^ 



1 — LU S H W 

T s = — — , H w = 

1 + to s n w 

where D w is the standard Wilson Dirac operator plus a negative parameter 
— m (0 < m < 2), and {u s } are a set of weights specified by an exact formula 
such that D c possesses the optimal chiral symmetry |2Zj- Since 

(D c + m,)- 1 = (1 - rmgy^D^img) - r], i 



2mo 
where 



D(m q ) = m q + (mo - m,/2)[l + 75 S (H w )], 

thus the quark propagator can be obtained by solving the system D(m q )Y = 1 
with nested conjugate gradient [22], which turns out to be highly efficient (in 
terms of the precision of chirality versus CPU time and memory storage) if 
the inner conjugate gradient loop is iterated with Neuberger's double pass 
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algorithm [3U] • For more details of our scheme of computing quark propagators, 
see Ref. [H]. 

We generate 100 gauge configurations with Wilson gauge action at — 6.1 
on the 20 3 x 40 lattice. Then we compute two sets of (point-to-point) quark 
propagators, for periodic and antiperiodic boundary conditions in the time 
direction respectively. Here the boundary condition in any spatial direction is 
always periodic. Now we use the averaged quark propagator to compute the 
time correlation function for any hadronic observable such that the effects due 
to finite T can be largely reduced [32] . 

Fixing m = 1.3, we project out 16 low-lying eigenmodes of \H W \ and 
perform the nested conjugate gradient in the complement of the vector space 
spanned by these eigenmodes. For N s = 128, the weights {uj s } are fixed with 
X min = 0.18 and X max = 6.3, where X min < X(\H W \) < X max for all gauge 
configurations. 

For each configuration, (point-to-point) quark propagators are computed 
for 30 bare quark masses in the range 0.03 < m q a < 0.8, with stopping criteria 
10~ n and 2 x 10~ 12 for the outer and inner conjugate gradient loops respec- 
tively. Then the chiral symmetry breaking due to finite N s (= 128) is less than 

io- 14 , 



a 



Y \S 2 Y _ 



yty 



< 10 



-ii 



for every iteration of the nested conjugate gradient, and the norm of the resid- 
ual vector for each column of the quark propagator is less than 2 x 10 -11 



\\(D c + m q )Y - 1\\ < 2 x 10 



-ii 



3 Determination of a and m s 

After the quark propagators have been computed, we first measure the pion 
propagator and its time correlation function, and extract the pion mass (m^a) 
and the pion decay constant (f n a). With the experimental input f w = 132 
MeV, we determine a" 1 = 2.237(76) GeV. 

The bare mass of strange quark is determined by extracting the mass of 
vector meson from the time correlation function 

1 3 

C v(t) = EE^lft + m q ) x ^{D c + m q ) l) 

6 fM=l X 

At m q a = 0.08, M v a = 0.4601(44), which gives M v = 1029(10) MeV, in good 
agreement with the mass of 0(1020). Thus we take the strange quark bare 
mass to be m s a = 0.08. Then we have 10 quark masses smaller than m s , 
i.e., m u a = 0.03,0.035,0.04,0.045,0.05,0.055,0.06,0.065,0.07,0.075. In this 
paper, we work in the isospin limit m u = m^. 
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5 10 15 20 2 4 6 8 10 12 14 16 

t t 

(a) (b) 

Figure 1: (color online) (a) The eigenvalue A + (t) of the lowest positive parity 
state, for m u a = 0.03. The solid line is the single exponential fit for 9 < t < 13. 
(b) The effective mass M eff (t) = \n[A(t)/A(t + 1)] of A+(t) in Fig. Ek, 



4 The 3x3 correlation matrix for Q(ududs) 

Next we compute the propagators ((O \) xa {0 y g) with fixed y = (0,0), and 
their time correlation functions Cfj(t) with ± parity 

C±(t) = fetr [i^(O i (f,t)O i (0,0)> / 

where the trace sums over the Dirac space, and the subscripts / and U denote 
fermionic average and gauge field ensemble average respectively. Then the 
3x3 correlation matrix C ± (t) = {C^(t)} can be constructed. Now, with 
the variational parameter to, we diagonalize the normalized correlation matrix 
C ± (^o)~ 1/2 C ,± (^)C ,± (^o)" 1/2 and obtain eigenvalues {Af (*)}, and from which to 
extract the masses {mf } of the lowest-lying and two excited states for ± parity 
respectively. In general, by varying to, one could minimize the errors of the 
masses extracted from the eigenvalues, as well as to disentangle (optimally) the 
lowest-lying states from the excited ones, as shown in Refs. ^21 E3] • However, 
in this case, the relevant quantities (e.g., the effective masses) extracted from 
unnormalized correlation matrix seem to be as good as those of the normalized 
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(a) (b) 

Figure 2: (color online) (a) The eigenvalue A~(t) of the lowest negative parity 
state, for m u a = 0.03. The solid line is the single exponential fit for 11 < t < 
17. (b) The effective mass of A~(t) in Fig. |2^. 



ones. Thus we restrict to unnormalized C ± (t) in the following. Then the mass 
mf can be extracted by single exponential fit to Af(t), for the range of t in 
which the effective mass M e ff(t) = \n[A(t)/A(t + 1)] attains a plateau. 

In Figs, n - IH the eigenvalues A ± (t) corresponding to the lowest-lying 
and first excited states with J p = l/2 ± are plotted versus the time slices, 
for m u a = 0.03 (the smallest quark mass in this study), together with their 
effective mass plots. Here we have suppressed any data point which has error 
(jackknife with single elimination) larger than its mean value. In each case, 
the mass m can be extracted by single exponential fit to A ± (t), for the range 
of t in which the effective mass M e ff(t) = \n[A(t) / A(t + 1)] attains a plateau. 
The results are (in units of a -1 ): 

™(l/2 + )i owes t_l ying = 1-034(80) 

^(V2 + )i st excited = 1-505(137) 
m(l/2-) lowegt . lying = 0.8045(23) 

ro(V2-)lst excited = 1-190(388) 

where all fits have confidence level greater than 0.6 and x 2 /d-o.f. < 1. Ob- 
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(a) (b) 

Figure 3: (color online) (a) The eigenvalue A + {t) of the first excited state with 
J p = l/2 + , for m u a = 0.03. The solid line is the single exponential fit for 
7 < t < 11. (b) The effective mass of A + (t) in Fig. Hi. 



viously, the lowest-lying and the first excited states are disentangled in both 
parity channels. 

In Fig. El the masses of the lowest-lying states with J p = 1/2^ are plotted 
versus m%. An interesting feature emerges in the positive parity channel. 
Its mass starts to fall more rapidly around the regime m u a ~ 0.045 (i.e., 
m u ~ 0.56m s ), signaling the onset of certain attractive interactions which 
lower the energy of this pentaquark state. Even taking into account of the 
error bar, the signal is unambiguous, as shown in Fig. \I\ior m u a = 0.03 (the 
smallest quark mass in this study). We suspect that this is the manifestation 
of diquark correlations when m u ^ becomes sufficiently small. To check, we 
compare it with the positive parity states extracted from the time-correlation 
functions of 0\, O2, and O3 respectively (see Figs. I9TITT1 in the Appendix). 
Then we see that similar phenomenon also happens in the positive parity 
channel of 3 (Fig. [Uj), but not in O x (Fig. EJ) or 2 (Fig. HUJ). Since 3 is the 
diquark-diquark-antiquark operator, it is consistent with our interpretation 
that diquark correlations emerge when m u a < 0.045 ~ 0.56m s . On the other 
hand, one may wonder whether it could be an artifact due to low statistics. 
Our argument is that if it were due to low statistics, it must also appear in the 
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(a) (b) 

Figure 4: (color online) (a) The eigenvalue A~(t) of the first excited state with 
J p = 1/2", for m u a = 0.03. The solid line is the single exponential fit for 
9 < t < 13. Note that the error becomes very large for t > 13. (b) The 
effective mass of A~ (t) in Fig. Ex- 



positive parity channel of 0\ and O2. However, it is not the case. So we rule 
out the possibility that this rapid decrease of the mass of the positive parity 
state in the regime m u a < 0.045 ~ 0.56m s is due to low statistics. In other 
words, it suggests that 3 has the largest overlap with the pentaquark state, 
and the diquark correlations may play an important role in forming @ + . 

The next question is how to perform the chiral extrapolation to the physi- 
cal limit where m n = 135 MeV (m u ^ ~ m s /25). From the viewpoint of chiral 
perturbation theory, one should use the set of data points with the smallest 
m u (m^). Moreover, since our m™ n is about 440 MeV, which may not be suf- 
ficently small to capture the chiral log behavior in chiral perturbation theroy, 
thus we only use the lowest order terms (i.e., linear in m\) for chiral extrapola- 
tion. Observing the onset of diquark correlations around m u a ~ 0.45, we nat- 
urally pick the smallest four masses (i.e., with m u a = 0.03, 0.035, 0.04, 0.045 ~ 
0.56m s ) for chiral extrapolation linear in m^. At m n = 135 MeV, we obtain the 
masses of the lowest-lying states: m(l/2~) = 1424(57) MeV, and m(l/2 + ) = 
1562(121) MeV, which agree with the masses of rrtK + ^jv — 1430 MeV, and 
0(1540) respectively. For the positive parity state, we have also performed 
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a fully correlated fit with the smallest four masses, employing the procedure 
adopted in Ref. Our result is 1554(150) MeV with x%u/ d -°-f- = 1-17(24), 
in agreement with the result of uncorrelated fit. By varying the fitting range 
of t, and the number of mass points, we estimate the systematic error to be 
180 MeV. 
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Figure 5: (color online) The masses of the lowest lying states of Q(ududs), 
extracted from the eigenvalues of the 3x3 correlation matrix of 0±, O2 and 
O3. The solid lines are chiral extrapolation (linear in m\) using the smallest 
four masses. 



5 Distinguishing the KN scattering states 

Now the question is whether the lowest-lying J p = l/2 ± states extracted 
from the 3x3 correlation matrix are scattering states or resonances. In order 
to obtain the mass spectrum of KN scattering states (which excludes any 
pentaquark states), we consider the time-correlation function of KN operator 
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Figure 6: (color online) The masses of the lowest-lying KN scattering states 
extracted from CxNif) ©? versus the naive estimates (with dotted lines). The 
solid lines (for J p = 1/2") are chiral extrapolations using the smallest four 
masses. 



without any exchange of quarks between K and iV in its propagator, i.e., 
the interaction between K and N is only through the exchange of gluons. 
Explicitly, 



C 



1 T74 



(N(x,t)N(0,0)) f {K(x,t)K(0,0)) 



(5) 



where N = [u T C7sd]d, and K = S7511. Note that the parity projector has 
taken into account of the intrinsic parity of kaon. 

The masses of lowest-lying KN scattering states are plotted in Fig. |Hl 
versus the naive estimates. Here the naive estimates are + for the 
and Jm? K + (2%/L) 2 + Jm 2 N + (2ir/L) 2 for the p-wave, where L is the 



s-wave, 



lattice size in spatial directions, and m K and m N are masses extracted from 
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Figure 7: (color online) The mass of the lowest-lying J p = l/2 + state extracted 
from the 3x3 correlation matrix, versus the s-wave of KNrj' ghost state (with 
dotted lines). 



kaon and nucleon time-correlation functions respectively. 

For the J p = 1/2" state, using the smallest four masses for chiral ex- 
trapolation to = 135 MeV, we obtain m KN {l/2~) = 1433(72) MeV, in 
agreement with rrix + ttijv — 1430 MeV. Further, its mass spectrum is almost 
identical to that of the lowest J p = 1/2" state of G in Fig. El for the entire 
range of m u . Thus we identify the lowest J p = 1/2" state of Q(ududs) with 
the KN s-wave scattering state. 

On the other hand, for the J p = l/2 + state in Fig. El its mass is higher 
than the naive estimate \Jm 2 K + (2n/L) 2 + \Jm 2 N + (2tt/L) 2 . This suggests 
that the KN p-wave (in the quenched approximation) in a finite torus is more 
complicated than just two free particles with momenta px = ~Pn = ^TrSt/L, 
i.e., their interaction through gluon exchanges cannot be ignored. Further, the 
mass of KN p-wave scattering state in Fig. |H1 is always higher than the mass 
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Figure 8: (color online) The mass of the lowest-lying J p = l/2 + state extracted 
from the 3x3 correlation matrix, versus the KN* s-wave scattering state with 
mass m K + m N *. 



of the lowest-lying J p = l/2 + state in Fig. E| In particular, for m u < m s , the 
former is significantly higher than the latter. This seems to suggest that the 
lowest J p = l/2 + state of Q(ududs) is different from the KN p-wave scattering 
state. 

However, there are other two hadron scattering states which also have 
quantum numbers J p = l/2 + and S = +1, namely, the s-wave scattering 
state of KNrj' (where rj' is an artifact due to quenched approximation), and 
the s-wave scattering state of KN* (where N* is the negative parity state of 
nucleon). In the following, we check whether the lowest-lying J p = l/2 + state 
extracted from the 3x3 correlation matrix could possibily be any one of these 
two hadron scattering states. 

The s-wave of KNrj' has J p = l/2 + , and its mass is estimated to be 
rriK + rriN + [Tl]. In Fig. [7[ it is clear that the lowest-lying J p = l/2 + 
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state of is different from the s-wave of KNrj' ghost state. Otherwise, they 
should almost identical for the entire range of m u , as in the case of KN s- 
wave and the lowest-lying state with J p = l/2~ in Fig. El So we exclude 
the possibility that the lowest-lying J p = l/2 + state extracted from the 3x3 
correlation matrix is due to quenched artifacts. Note that Mathur et. al. 
has claimed that the quenched artifact KNrj' (with negative spectral weight) 
would be seen in the positive parity channel if m n < 300 MeV, for the 0\ 
[K ® 75iV) operator. If this claim holds for any cases (gauge actions, lattice 
fermions, operators, etc), then we should not see KNrj', since our smallest 
pion mass is larger than 400 MeV. 

Next, we turn to the s-wave of KN*, which also has J p = l/2 + , and 
its mass is estimated to be m^ + m^*. In Fig. |H1 it is clear that the mass 
of KN* s-wave is different from the mass of the lowest J p = l/2 + state 
extracted from the 3x3 correlation matrix, for the entire range of m u (except 
at the crossover m u a ~ 0.03). In the physical pion limit, the mass of the 
lowest J p = 1/2+ state is 1562(121) MeV, while m K + m N * is 2137(142) MeV. 
This rules out the possibility that the lowest-lying J p = l/2 + state extracted 
from the 3x3 correlation matrix is the KN* s-wave scattering state. Note 
that this disagrees with the claim of Ref. (221, which was based on their 
result obtained with Wilson quarks at unphysical quark masses. However, 
there is no reason to expect that baryon masses computed with two different 
lattice (fermion and gauge) actions would agree with each other, except in the 
continuum (and infinite volume) limit with physical quark masses. Thus, it 
is not surprising to see that, at unphysically large quark masses, our result of 
the mass of the lowest-lying state with even parity, m(l/2 + ), is higher than 
rax + itin*, while the result of Ref. [22] seems to suggest that m(l/2 + ) ~ 
mK + mN*. Now the emerging problem for Wilson fermion is to check whether 
the relationship m(l/2 + ) ~ mK + m^* remains valid in the physical pion limit. 
We suspect that even for Wilson fermion, it would also exhibit the inequality 
m(l/2 + ) < m^ + ?tijv* at sufficiently small quark masses, though the location 
of the crossover might be different from that of Fig. |H1 At this point, it is 
instructive to check the result of Mathur et al. [Tl|, which was obtained with 
the Overlap fermion. From Ref. fl] and hep-ph/0306199, one can extract 
m(l/2 + ), m^i and m^* respectively, then one immediately sees that m(l/2 + ) 
is larger than m^ + m^* for m\ > 0.27 GeV 2 , but becomes smaller than 
m^ + mjq* for m\ < 0.27 GeV 2 , and it tends to ~ 1650 MeV in the physical 
pion limit. Obviously, the result of Mathur et al. [H] also suggests that the 
lowest state with even parity is different from the KN* s-wave, in agreement 
with our conclusion. Now a (physically irrelevant) question is where the KN* 
s-wave scattering state lies, say, at unphysically heavy quark masses. Our 
conjecture is that in a finite torus, the KN* s-wave might turn out to be much 
heavier than the naive estimate mx + rriN* (similar to the case of KN p-wave 
which is much heavier than its naive estimate, as shown in Fig. thus it 
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always lies above the lowest J p = l/2 + state, as one of the excited states in 
the positive parity channel. 

Now, after ruling out the possibilities of being K N p-wave, KNrj' s-wave, 
or KN* s-wave, the lowest-lying J p = l/2 + state extracted from the 3x3 
correlation matrix seems to be nothing but a resonance. If it is identified with 
6 + (1540), then it predicts that the parity of + (154O) is positive. 

6 Discussions and concluding remarks 

Now we return to Table [T] to discuss what causes the different claims in these 
exploratory lattice studies. 

Now, back to Fig. |S| If we had not measured any data points with 
m u < m s /2, then we could not have seen the rapid decrease of mass of the 
positive parity state (in the regime m u a < 0.045), of which we interpret as 
the manifestation of diquark correlations at sufficently small m u ^- Conse- 
quently, chiral extrapolations using data points with m u > m s /2 must yield 
much higher masses, especially for the positive parity state. For example, 
if we use the data points with m u a = 0.05, 0.055, 0.06, 0.065 < m s a, then 
we obtain m(l/2~) = 1666(47) MeV, and m(l/2+) = 2178(104) MeV, with 
chiral extrapolation linear in m^. Further, if we use the data points with 
m u a = 0.07, 0.075, 0.08, 0.085(> m s a), then we obtain m(l/2~) = 1725(52) 
MeV, and m(l/2 + ) = 2616(103) MeV. In either one of these two cases, it 
seemingly rules out the possibility that the positive parity channel could ac- 
commodate any state with mass 1540 MeV. However, it is only an artifact due 
to chiral extrapolation with data points too far away from the physical reality 
(m U) d ~ m s /25). This explains why the masses of the positive parity state in 
Refs. [121121] are so high comparing to our result. Note that, in Refs. [T§1 |2"T] . 
the total number of data points is four, in which the number of "physical" 
data points (i.e. with m u < m s ) is only one and zero (see Table [TJ. Even 
though they claimed that the positive parity channel could not accommodate 
G + (1540), it is most likely just an artifact due to their chiral extrapolations 
with data points too far away from the physical reality. 

Next, we turn to the claims of Refs. fHl EJ- From our mass spectra of 0\ 
(Fig. EI) and O2 (Fig. ITU)) , one can explain why Csikor et al. ^H] (with 0\ and 
O2) obtained a rather high mass for the positive parity state, while Mathur 
et al. PI] (with 0\) did not see pentaquark resonance in the positive parity 
channel. Namely, the interpolating operators 0\ and O2 have little overalp 
with the pentaquark state in the positive parity channel. In fact, our mass 
spectra of 0\ and O2 are consistent with those of Refs. [TTfl I14j . 

Finally we discuss the results of this paper. This is the first lattice QCD 
study on G with 3x3 correlation matrix. Presumably, it should provide a more 
reliable answer to the questions of signal/parity of + than other lattice studies 
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with only one operator. However, this is a quenched lattice QCD calculation 
(like other lattice studies on + so far), with only one volume, and one lattice 
spacing, thus it is difficult for us to estimate the systematic error. For the 
lowest J p = l/2~ state, it is identified with the KN s-wave, by comparing 
its mass with + tuk- For the lowest J p = l/2 + state, by comparing its 
mass (as a function of m 2 ) to those of two-hadron scattering states having the 
same quantum numbers, it seems unlikely to be identified with any one of the 
following two-hadron scattering states: KN p-wave, KNt]' s-wave, and KN* 
s-wave. However, before it can be confirmed to be a resonance, it is necessary 
to check whether its mass and spectral weight are volume independent. To this 
end, we are performing computations on the 24 3 x 48 lattice with the same 
lattice spacing (i.e., Wilson gauge action at (3 = 6.1). If the lowest-lying state 
with J p = l/2 + turns out to be a scattering state, then there is no evidence 
of pentaquark resonance in our study. On the other hand, if it turns out to 
be a resonance, then it can be identifed with + (154O) since its mass is close 
to 1540 MeV. Nevertheless, one still has to find out whether its decay width 
could be as small as 15 MeV (compatible to that of + ), which is the most 
challenging problem pertaining to + . 

Appendix 

In this appendix, we present our results of the masses of the even and odd 
parity states extracted from the time-correlation functions of operators 0\ ((TJ), 
2 (0), and 3 respectively. 

In Figs. lUlim the masses of the J = 1/2^ states are plotted versus m%, for 
0\, 2 , and 3 respectively. Here all mass fits have confidence level greater 
than 0.6 and x 2 /d.o.f. < 1. For 0\ and O2, their masses vary smoothly with 
respect to m^., in both parity channels. On the other hand, the positive parity 
state of O3 undergoes a rapid decrease for m u a < 0.045, similar to the behavior 
of the lowest-lying J p = l/2 + state extracted from the 3x3 correlation matrix 
(see Fig. [5J). This seems to signal an onset of certain attractive interactions 
when the quark mass m u ^ becomes sufficiently small. We conjecture that this 
is the manifestation of diquark correlations when m u ^ approaches the physical 
limit. 

Following the same argument in Section 4, we pick the smallest four masses 
(i.e., with m u a = 0.03, 0.035, 0.04, 0.045) for chiral extrapolation (linear in m 2 ). 
At physical pion mass = 135 MeV, we obtain 





m(l/2~) 


= 1430(66)MeV, 


m(l/2+ 


) = 2301(134)MeV 


2 : 


m(l/2~) 


= 1430(67)MeV, 


m(l/2+ 


) = 2346(156)MeV 


O3: 


m(l/2 _ ) 


= 1446(71)MeV, 


m(l/2 + 


) = 1843(136)MeV 



Note that for the J p = 1/2 state, all three operators give almost the 
same mass which coincides with that [1424(57) MeV] extracted from the 3x3 
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correlation matrix of 0\, O2, and O3. However, for the J p = l/2 + state, the 
operator O3 (the diquark-diquark-antiquark operator motivated by the Jaffe- 
Wilczek model) gives the lowest mass, which is the closest to that [1562(121) 
MeV] extracted from the 3x3 correlation matrix of 0±, O2, and O3. This 
seems to imply that among {Oi,i = 1,2,3}, O3 has the largest overlap with 
the pentaquark state, and the diquark correlations may play an important role 
in forming O + (1540). 

Further, we also observe that the diagonalization of the 3x3 correlation 
matrix indeed disentangles the contributions of the excited states, and gives a 
smaller mass for the lowest-lying J p = l/2 + state. 
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Figure 9: (color online) The masses of the even and odd parity states ex- 
tracted from the time-correlation function of 0\. The solid lines are chiral 
extrapolation (linear in m^) using the smallest four masses. 
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Figure 10: (color online) The masses of the even and odd parity states ex- 
tracted from the time-correlation function of O2. The solid lines are chiral 
extrapolation (linear in m^) using the smallest four masses. 
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